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. . . Abstract 
CN ■ 

■ The (3+1) dimensional expansion of the quark gluon plasma (QGP) produced at finite baryon 
, density has been studied using relativistic hydro dynamical approach. The pressure functional of 

the equation of state (EoS) has been determined for the interacting nuclear matter with mesons 

■ exchange. The EoS has been used to solve hydrodynamical equations using RHLLE algorithm. 
The space time expansion of the plasma has been studied for the cases of SPS energy and RHIC 
energy both at finite baryon density and for a hypothetical case of SPS energy at zero baryon 
density. The space-time evolution is slowed and the life times of QGP and mixed phases are 
shortened in the presence of finite baryon density. The space time integrated total photon yields 

i have been estimated by convoluting the static emission rates with the space time expansion of 

I the plasma. It has been shown that the total photon yield at zero rapidity is not significantly 

lO ■ affected by the baryon density for SPS energy. The total photon yield is unaffected by the 

, Landau-Pomeranchuk-Migdal effect at SPS energy at zero baryon density, as the quark matter 

' contribution to photon yield is less compared to hadronic matter. 

(N ' 
O 

QJ , Electromagnetic processes such as photons and dileptons production are important to identify 

the quark gluon plasma (QGP) expected to be formed in the ultra relativistic heavy ion collisions 
at SPS, RHIC and LHC energies ^, §]. Consequent to the availability of the experimental 
^ . data reported by WA98 collaboration ||5[ of single photons from Pb+Pb collisions at CERN SPS, 

^ I these studies gained much interest in recent times. The experimental spectra however include 

photons coming from various sources like decays, prompt and thermal photons. Photons are 
produced at each stage from an expanding volume of plasma being both in QGP and hadron 
phases. The yields of these electromagnetic signals are quite sensitive to the space time evolution 
of the QGP. Assuming a first order phase transition, the plasma will expand and evolve through a 
QGP, mixed and hadron phases governed by the appropriate equation of state. Thermal photons 
are emitted at a later stage which depends on the space time evolution scenario considered 
consequently, the expansion dynamics determines the observed spectra. Initially, the plasma 
expands longitudinally, however at high energy densities and when plasma life time is large, 
the transverse expansion effects become significant and can not be ignored |^]. Further, the 
experimental data indicates the presence of finite baryon density at freezeout temperature. The 
evolution becomes more complicated if the plasma is unsaturated or formed at finite baryon 
density. These physical conditions of the plasma affect both the basic photon production rate 
and also the total space time integrated photon yields. In the present work we incorporate the 
finite baryon density effects in both the production rates and the plasma expansion dynamics. 

The photon production rates from QGP phase have been investigated in detail in These 
studies include the production rates at effective one loop and two loop levels. Further, multiple 
soft scatterings of the fermion during photon emission reduce the photon coherence lengths. 



known as Landau-Pomeranchuk-Migdal (LPM) effect, and suppress the production rates |12 
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15, |T6|. The photon rates from the hadron phase were estimated considering the comprehensive 



set of processes in 18], and the prompt decays were estimated in |19, 2C]. Several groups 
have already reported the space-time integrated photon yield calculations for SPS, RHIC and 
LHC energies considering (i) various sets of equations of state, (ii) considering one or three 
dimensional space time expansion of plasma, (iii) with or without boost invariance, for details 
see |,|,|,|,| 0,11. 

In the following, the expansion dynamics of the plasma is determined by the EoS and the rela- 
tivistic hydrodynamical equations. The conservation equations for the energy {E), momentum 
(M) and net baryon number (R) densities are given by, 

di.T'"' = and, df.N^' = (1) 

E = T^o = (e+p)j''-p (2) 

M = = (e + p)-f^v (3) 

R = N'^ =nu° = n-f (4) 

Here, = 7(1, v) is the fluid velocity. In the above, e,p,n are the energy, pressure and net 
baryon number densities in the local rest frame of the fluid. The conservation equations for 
E, M, R of Eqs. (2-4) for cylindrical geometry with the longitudinal (z-axis) boost invariance 
and at z=0 are as discussed in ||2^. The shock just at the QGP formation time is represented by 
the Riemann initial value problem. The full solution to the Riemann problem is approximated by 
the Godunov scheme consisting of a rarefaction wave and the contact discontinuity. This scheme 
approximates the Riemann shock propagation by assuming a region of constant flow between 
the discontinuous surfaces. Our present calculations have been performed using Relativistic 
HLLE (RHLLE) program to solve the hydrodynamical equations |2^. An important input 
for this method is the equation of state (EoS) specifying the pressure as a functional of the 
energy density and the number density p{e,n). Once an EoS is specified, the hydrodynamical 
equations uniquely determine the expansion dynamics. We obtain the EoS for hadronic, mixed 
and QGP states at finite baryon density. In this study, the EoS for the QGP state is taken from 
bag model consisting of u,d quarks and gluons. The pressure, energy and entropy densities are 
given by, 

^ 377r2 1 - B (5) 

dP 

p= {e- 4B)/3 ; s = — ; e = Ts + fin-p (6) 

In above, the fi is the baryo chemical potential representing the finite baryon density effects. The 
hadronic EoS is obtained from a general class of thermodynamically self consistent equations of 
the interacting nuclear matter that reproduce the ground state matter properties. The present 
EoS includes the a,uj mesons interacting with nucleons and anti nucleons (i.e. masses upto 1.0 
GeV) and a free gas of vr, p, r/ mesons. The hadron pressure (ph), net baryon number density 



(n), the scalar density (p) and the partial pressures of baryon and mesons {pN,Pi) are as in |25| 
given by, 

rn rps 

Ph{T, ^)=pN + EiPi + ni/(n) - / i^{n')dn' - p,S(p,) + / S(p'Jdp', (7) 

Jo Jo 

n{T,^x) = ^J d'k[fT{Elv)-fT{El-u)], with u = ^Ji-v{n) (8) 
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(2^3) 



M* 
'Et 



[fT{Elv) + fT{El-v)] 



El = Jk^ + M*^; M* = M-S{ps) ; iy{n) = C^n - Cjps ; S{ps 



PN =T 



9n 



(27r3) 
Pi{T;mi) 



d^k 



-T 



ln(l + e-(^^*-^)/^) + ln(l + g-^^^+^^Z^ 
ln(l - e-^V^/T)) 
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(27r3) 



d-'q 



(9) 

(10) 
(11) 

(12) 



In the above, fx is the fermi distribution at temperature T and baryo chemical potential p 
representing the baryon density, rrii are the meson masses, M is the nucleon bare mass and 
gi,gN are the meson and nucleon spin-isospin degeneracies. The energy density (e) and entropy 
(s) are obtained using thermodynamic relations as in Eq. 6. The effective chemical potential 
i'(n) and the effective mass M* are obtained using the parameters given in = 238.08Gey~^, 
C2 = 296mGeV-^, Cj = 0.183 |5|. The Eqs.(8,9,10) are solved seff consistently . The phase 
boundary in T, /i plane, shown in Fig. 1(a), is obtained from the Gibbs conditions given by, 



Ph = Pq, Th = Tg, flh = Pq (13) 

The mixed phase for a given (e,n) is obtained from, 

e = \qeq{T*) + {l-\)eh{T*) and n = \nq{T*) + {I - \q)nh{T*) (14) 

In the above, the subscripts h denote the qgp,hadron phases. \q is the fraction of volume 
occupied by the QGP in the mixed phase. The EoS specifies the pressure as a functional of 
p(e,n) in the required region in (e,n) plane as shown in the Fig. 1 (nearly same as Fig. 1 of 
||2^ ). Fig. 1(a) shows the transition temperature and chemical potential given by Eq. 13. The 
transition temperature at zero baryon density is 169.9 MeV. Fig. 1(b) shows the energy density 
variation and the Fig. 1(c) shows the number density variation along the curve of Fig. 1(a). 
Fig. 1(d) shows the energy density and number density for both quark and hadron phases for T, p 
values along the curve of Fig. 1(a). The hadron phase is the narrow strip in the semi- infinite 
(e,n) plane as shown in Fig. 1(d) by the mixed phase above and below by the T = line of 
hadron phase. However, owing to small life times of QGP and mixed phases (10-20 fm/c), these 
phases occupy very small region of space time evolution diagrams which will be discussed later 
(Fig. 2). The T* of Eq. 14 are determined using Figs. l(b,c,a). The constants used in Fig. 1 
are eo = O.U77n^GeV^ and no = O.Wh^GeV^. 

We studied the expansion dynamics at finite baryon density for two cases; case (i) for SPS energy 
and case (ii) for the RHIC energy. The case (iii) for SPS energy at zero baryon density is studied 
to compare the effects of finite baryon density. The initial energy density and the number density 

at time r = Tj have been assumed to be of wounded nucleon form, as given by /(r) = — 
and R is the nuclear radius taken to be 6.0 fm. The formation times are taken as Ti=l fm/c 
for SPS energy and rj=0.6 fm/c for the RHIC energies respectively. The parameters for RHIC 
energies were taken from [^3|. The entropy (s), energy and number densities corresponding to 
SPS case of 158 AGeV of Pb+Pb collision is obtained by solving the simultaneous equations 
given by. 



TtR UbTi 



dNB 
dy 

irK STi = 3.6—— 
dy 
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We have taken = 80, and ^ = 700 from the experimental data, which corresponds to a 
baryon density of n{,/no=4.5 . These space-time evolution of the plasma for the cases of SPS 
energy at finite baryon density is shown in Fig. 2 for different phases as indicated in the figure 
in various colors. The evolution of zero baryon density is also shown in Fig. 2, by thick colored 
lines representing the boundary of corresponding phases. The radial extension and the proper 
time (at z=0) are in units of nuclear radius. The outermost envelope of the hadron phase in 
Fig. 2 corresponds to the freeze out isotherm taken to be 120.0 MeV. The QGP phase region is 
almost the same for the two cases (slightly less at finite density). It can be seen in Fig. 2 that at 
SPS energy the baryon density decreases the life times of QGP and mixed phases of the plasma 
as well as the freeze out time. The transverse extension is also affected by the baryon density. 
Similar cases of space time evolution have been studied for RHIC energy. At RHIC energies, 
the life times are considerably larger together with the large transverse extension (see Fig. 2 of 
125 



The temperature, chemical potential, constant energy density and the transverse velocity pro- 
files in the space time, as the plasma expands, have been studied for these cases. The initial 
temperatures are larger in the QGP phase which in time go over to the constant temperature 
profiles corresponding to mixed phases. Subsequently, the temperature falls very rapidly as the 
plasma cools. The maximum temperatures (at r=0 for wounded nucleon form) reached are ap- 
proximately 239.06 MeV, 244.25 MeV for SPS energies for the cases of with and without baryon 
density respectively and 348.4 MeV at RHIC energy. The average energy densities correspond 
to the temperatures of 214 MeV, 219.5 MeV at SPS energy for the cases of with and without 



baryon density. These results are comparable to those of Table 1 of |21]. 



The energy fiow in space-time while in expansion is demonstrated by the constant energy con- 
tours. These contours for the SPS cases have been compared in Fig. 3. The solid curves are 
with baryon density and the dashed curves are for zero density. The constant energy density 
values are as mentioned in the figure in units if GeV/fm^. It can be seen from the figures that 
the energy flow gets slower in the presence of baryon density. Similarly, the constant transverse 
velocity proflles for various values of /? = v/c have been shown for these cases in Fig. 4. 

In the following we discuss the static rates for photon emission needed to calculate the total 
photon yields. The annihilation processes {qq gj) and the QCD Compton {qg — > 57, qg — > ^7) 
processes have been considered for the thermal photon production at effective one loop level 
1^, ^. The photon emission rates from quark matter at one loop have been taken from the 



Traxler et al. Q, and at two loop level for non zero chemical potentials as reported in |17| 
We reported the photon production rate upto two loop level for a general case of chemically 
unsaturated plasma at flnite baryon density in |l^. The bremsstrahlung and annihilation with 
scattering (aws) processes at two loop level have been considered. At finite baryon density, 
the photon production at one loop level is dominated by the compton process compared to 
annihilation process. At two loop level, the bremsstrahlung radiation is affected whereas the 
aws process is insensitive to the baryon density presence. The bremsstrahlung radiation from 
quark is enhanced and from anti quark is suppressed at finite baryo chemical potential. 

The photon emission rates from the hadron phase is determined by considering various meson 
reactions, such as p — > tttt^, vrp — > tt'j, vrvr — > p^, lo n^. The vrp tt^ estimate also includes 
the ai intermediate resonance. The numerical results of a detailed study of these decays 
has been fitted by analytical expressions as a function of temperature and photon energy, as 
reported in ||2^, 27 , [l| . In the present work, we used these analytical formulae for the photon 



rate calculations. It has been shown that from temperature in the range of 100-200 MeV, these 
formulae reproduce quite well the numerical results from detailed calculations for the photon 
energies from 200 MeV and upto 5.0 GeV. The one dimensional expansion of the plasma with 
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boost invariance has been already studied with these emission rates with and without baryon 
density. It was found that the finite baryon density has no effect on the total photon yields |28|. 
The total photon yield integrated over the plasma space-time hyper volume at zero rapidity is 
given by, 



(15) 

In above, the subscripts q and h denote the photon emission rates from QGP and hadron 
phases. The Tf and denote the freeze out (outermost) boundary shown in Fig. 2. The 
rates correspond to the photons emitted in a fluid rest frame. In order to compare with the 
experimental photon yields, the transverse velocity distributions of Fig. 4 have been used to 
transform the photon energy as given by E^ = p^u^. Fig. 5 shows the photon transverse 
momentum rapidity density distributions at zero rapidity for the 158 A GeV SPS energy for 
Pb+Pb collision compared with WA98 data. The finite and zero density cases are shown in 
figure. As seen in figure, the finite baryon density has very little effect on the photon spectrum, 
though the yield seems to be reduced marginally. This result is surprising considering that the 
expansion dynamics is affected as shown in Fig. 2. However, at these SPS energies, the overall 
plasma life time is small and the expansion effects and modified emission rates due to finite 
chemical potential cancel each other. However, the fit to experimental data is slightly improved 
when the finite density is considered (see Fig. 5). The prompt contribution in the present study 



has been taken from [^]. It has been noticed that the prompt estimates of |2^ and [19| differ by 
about a factor of two together with roughly same slope in the region of 1.5-4 GeV. Therefore, if 
one considers the prompt contribution twice as large as given in |19|, our calculated spectrum at 
finite density matches well with the experimental data. At RHIC energies, the baryon number 
densities are smaller and the finite density effects can be completely ignored. All these studies 
have been repeated for SPS case for different formation times rj=l.l fm/c,1.2 fm/c and the 
results are very much similar. The results for photon yields have not changed significantly with 
different formation times studied. The present study was partly reported in [p9|. 



The photon rates from QGP phase at two loop level of |17] are identical to the results of 
Aurenche et al. [11| for the equilibrated plasma, corrected by a factor of four. These rates and 
the total photon transverse momentum distributions shown in Fig. 5 for zero baryon density 
donot include the LPM effects. The effective two loop contributions to photon production rate 
from QGP phase consisting of bremsstrahlung and aws processes with and without LPM effect 
are discussed in |15, |l^, especially as shown in Fig. 7 of Ref. |16|. We have studied the 
consequences of photon rate inhibition due to LPM effect on the space time integrated photon 
yields for the case of zero density. For this purpose, we used the phenomenological expressions 
of |16| that include LPM effects. The ratio of total photon yield arising from the two loop 
processes with and without LPM effect been estimated and shown in Fig. 6. It has been found 
that the LPM effect has not changed the photon transverse momentum distribution, due to the 
fact that at SPS energy the quark matter contribution to photon yield is very less. The LPM 
effect becomes important at higher energies as the aws contribution to photon emission is more 
suppressed. The total quark matter contribution from one and two loop processes without LPM 
effects are calculated and the relative percentage contribution of one loop processes and the two 
loop processes have been obtained for the space-time integrated yields. It should be noted that 
the temperature and the chemical potential dependence has been integrated out through the 
space time integrations of Eq. 15 over the plasma volume. These are also shown in Fig. 6 where 
short dashes represent oneloop/total quark matter emission and long dashes similarly for the 
two loop contribution. The different color curves show the results for RHIC energy. At RHIC 
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energies the LPM effects are important as the QM contribution is important and therefore the 
affects the calculated yields. 



Conclusion 

The effects of finite baryon density and the (3+1) dimensional expansion of the quark gluon 
plasma has been studied for SPS and RHIC energies. The transverse expansion has been ob- 
tained following relativistic hydrodynamical approach and the relativistic HLLE method. The 
equation of state for pressure functional has been obtained considering various mesons and in- 
teracting nucleons. The space time evolution of the QGP is shown to be affected by the baryon 
density. The radial and the temporal extension of the plasma is shorter in the presence of finite 
baryon density. The total plasma space time integrated photon yields have been compared for 
various cases in the presence of baryon density. Calculations at SPS energies with and without 
baryon density showed that the total photon yield is not very sensitive to the baryon density, 
however individual contributions to photon yields from various phases differ. The LPM effect is 
also found to be not affecting the photon yields for SPS energies, where as at RHIC energy this 
effect is important. 
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Figure 2: The space-time expansion of the QGP for SPS case. The figure is for finite baryon 
density, and the corresponding colored sohd curves represent the extent of that phase for zero 
baryon density case. 
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Figure 3: constant energy contours for the SPS energy. The sohd curves are with density effects, 
and the dashed curves are for zero density case. 
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Figure 4: constant velocity contours for the SPS energy. The sohd curves are with density 
effects, and the dashed curves are for zero density case. 
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Figure 6: The suppression factor, for the two loop processes consisting bremsstrahlung and aws 
processes from QGP phase, arising from LPM effects. The LPM effect is shown by solid curves, 
the dotted curve shows the relative one loop contributions in total quark matter for the photon 
yields and similarly the dashes are for the relative two loop contribution. 
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